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We use a sample of 384 million BB decays collected with the BABAR detector at the PEP-II asymmetric e + e — storage 
ring to study the flavor-changing neutral current decays B — » K^t^i~ , where l + £~ is either e+e~ or /Lt+ fi~ . We 
present measurements in two dilepton mass bins, one below the J/tj} resonance and the other above, of the lepton 
forward-backward asymmetry Afb and the longitudinal K* polarization Fl in B — > K*l+l~ , along with isospin rate 
asymmetries in B — > K*£+£~ and B — * Kt^l~ final states. 

"66 : 
o ■ 

The decays B -> K^k + £~, where £ + £~ is either e+e" or , arise from flavor-changing neutral current processes 

that are forbidden at tree level in the Standard Model (SM). The lowest-order SM processes contributing to these 
O ■ decays are the photon penguin, Z penguin and W + W~ box diagrams. Their amplitudes are expressed in terms of 
hadronic form factors and effective Wilson coefficients CV, C9 and C10, representing the electromagnetic penguin 
£NJ diagram, and the vector part and the axial- vector part of the Z penguin and W + W~ box diagrams, respectively l[. 

New physics contributions may enter the penguin and box diagrams at the same order as the SM diagrams [2|. 
1— 1 We present measurements of the lepton forward-backward asymmetry Afb and longitudinal K* polarization Fl in 
B — > K*£ + £~ , along with isospin rate asymmetries in B — > K*£ + £~ and B — > K£ + £~ final states, in two bins of 
dilepton mass squared q 2 = m 2 +g _ , one below the J/ip resonance and the other above. Sensitive indirect searches for 
new physics effects using these observables Q are possible as hadronic uncertainties in calculations are expected to 
■ cancel [J]. BABAR results on branching fractions, direct CP and lepton-flavor asymmetries, Afb and Fjj have been 
previously published, however, only a limit for Afb in the low q 2 region was established [5(. 
CN ■ The CP-averaged isospin asymmetry 

: kW _ B{B° -> K^°£+£-) - rB(B ± -> K^^^) 

I 1 ~ B(B° -» K(*)°£+£-) + rB(B ± -> K(*)±£+£~) ^ 

where r = tq/t+ = 1/(1.07 ± 0.01) is the ratio of the B° and B + lifetimes Q, has a SM expectation of +6 — 13% 

as q 2 — * for B — > K*£ + £~ [3|. This is consistent with the measured asymmetry of 3±3% in £? — > K*j [J. ^A 

if* _ _n nn+0.005 



calculation of the predicted K* + and SM rates integrated over the low q 2 region gives Af* = — O.Q0^ ; u06 
Given that the expected SM isospin asymmetry arises from pure photon penguin contributions to B — > K*£ + £~, 
^ . there is no expectation of such an asymmetry in K£ + £~. In the high q 2 region, although there may be possible 

contributions from higher charmonium resonances, the measured asymmetry in J/ipK 1 -*" 1 is only a few percent 
;_i ' and any SM Af can be reasonably expected to be similarly insignificant. We measure isospin asymmetries in a 
low q 2 region 0.1 < q 2 < 7.02GeV 2 /c 4 , and a high region q 2 > 10.24 GeV 2 /c 4 , where any likely contributions from 
J/ip and ip(2S) have been removed by vetoing events 7.02 < q 2 < 10.24 GeV 2 /c 4 and 12.96 < q 2 < 14.06 GeV 2 /c 4 , 
respectively. 

The K* longitudinal polarization fraction Fl can be determined from the distribution of the angle cos 9k between 
the K and the B directions in the K* rest frame using a fit to cos 9k of the form [l^] 

^F L cos 2 9 K + - A {l-F L )(\-cos 2 9 K ) (2) 

Likewise, the lepton forward-backward asymmetry Afb can be determined from the distribution of the angle cosf^ 
between the £ + {£~) and the B(B) direction in the £ + £~ rest frame using a fit to cos9( of the form [l3| 

-F L {1 - cos 2 Bt) + |(1 - F L ){1 + cos 2 B t ) + A FB cos 9 t (3) 
4 

These angular measurements are also done in low and high q 2 regions, with the low region slightly more narrowly 
defined than above, 0.1 < q 2 < 6.25GeV 2 /c 4 , in order to absolutely remove any possible contributions from J/ip and 
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ip(2S) decays, which have a distinct angular structure but are otherwise indistinguishable from signal decays. As 
shown in Eq. 3, an experimental determination of F L is required to obtain a model- independent Afb result. 

Variations in Fl and Afb as a function of q 2 result from interference among the different amplitudes. The expected 
SM behavior of Fl and Afb, along with variations due to non-SM opposite-sign Wilson coefficients, is shown by 
the curves in Fig. JT] At low q 2 , where CV dominates, Afb is expected to be small with a zero-crossing point at 

1 constraint on the magnitude of CV comes from the measured inclusive 



q 2 ~ 4 GeV 2 /c 4 
branching fraction for b — > 

q 2 in the SM, the product C9C10 is expected to give a large positive asymmetry. Right-handed weak currents would 
have an opposite-sign C9C10, leading to a negative Afb contribution at high q 2 . Contributions from other non-SM 
processes could change the magnitudes and relative signs of C7, Cg and C10, and introduce complex phases between 
them [HI El I . 
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which corresponds to the limit q 2 — > 0, but its sign is unknown. At high 
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Figure 1: (a) Afb and (b) Fl results showing comparisons with SM (solid), CV = — C7 (long dash), C9C10 = -C7| M Cf M 
(short dash), C7 = —Cf M ,CgCio = — Cjf M Cf M (dash-dot). Statistical and systematic errors are added in quadrature. 
Expected Fl values integrated over each q region are also shown. ™- n - r '-- - nSMnSM 

identical to the two curves shown. 



The F L curves with C 9 C W = -C$ M C? M are nearly 



We use a data sample of 384 million BB pairs collected at the T(AS) resonance with the BaBAR detector [3] at the 
PEP-II asymmetric-energy e + e~ collider at SLAC. Tracking is provided by a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) in a 1.5-T magnetic field. We identify electrons with a CsI(Tl) electromagnetic 
calorimeter, and muons using an instrumented magnetic flux return. Electrons (muons) are required to have momenta 
p > 0.3(0.7) GeV/c in the laboratory frame. We combine photons with electrons when they are consistent with 
bremsstrahlung, and do not use electrons that are associated with a photon converting to a low-mass e + e~ pair. We 
identify K + using a detector of internally reflected Cherenkov light, as well as ionization energy loss measurements 
from the DCH and SVT. Charged tracks other than identified e, p and K candidates are treated as pions. Neutral 
pion candidates are formed from two photons with laboratory energies Ey > 50 MeV and an invariant mass between 
115 and 155MeV/c 2 . 

We reconstruct signal events in ten separate final states containing an e + e~ or pair, and a Kg(—* n + Tr~), 

K + , or K*(892) candidate with an invariant mass 0.82 < M(Ktt) < 0.97GeV/c 2 . We reconstruct K*° candidates 
in the final state K + ir~, and K* + candidates in the final states K + ir° and K®tt + (charge conjugation is implied 
throughout). Neutral K® — > ir + ir~ candidates are required to have an invariant mass consistent with the nominal 
K° mass [17j ] , and a flight distance from the primary interaction point which is more than three times its uncertainty. 
We also study final states K^h^/j^ , where h is a track with no particle identification (PID) requirement applied, to 
characterize backgrounds from hadrons misidentified as muons. Signal decays are characterized using the kinematic 
variables tties = a/s/4 — p*^ and AE = E* B — \i r s/2, where p* B and E B are the B momentum and energy in the 
T(4S) center-of-mass (CM) frame, and y^s is the total CM energy. We define a fit region tties > 5.2GeV/c 2 , 



with -0.07 < AE < 0.04 (-0.04 < AE < 0.04) GeV for e 
-0.08 <AE< 0.05 (-0.05 < AE < 0.05) GeV for high q 2 . 



(p + p ) final states in the low q 2 region, and 
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The main backgrounds arise from random combinations of leptons from semileptonic B and D decays, which are 
suppressed using event shape variables, vertexing information and missing energy combined in event selection neural 
networks (NNs). We use simulated samples of signal and background events in the construction of the NNs and, 
assuming rates consistent with accepted values |6| , we optimize the NN selections for best statistical significance in the 
number of expected signal events in our dataset. A further background contribution comes from B — » D(—> K^*'tt)tt 
decays where both pions are misidentified as leptons. This background is significant only in dimuon final states, 
as the pion misidcntification rate for electrons is < 0.1%. We veto these events by assigning the pion mass to a 
muon candidate and requiring that the invariant mass of the hypothetical K^n system be outside the range 1.84- 
1.90GeV/c 2 . We perform blind fits for both the isospin and angular observables, in which all event selection criteria 
were determined prior to examining events in the fit region. 

We directly fit the data with Af ( 1 as a floating parameter using a simultaneous unbinned maximum likelihood 
toes fit across all modes contributing to a particular Aj measurement. An ARGUS shape [3] with floating shape 
parameter and normalization is used to describe combinatorial background. For the signal, we use a fixed Gaussian 
shape unique to each final state, with mean and width determined from fits to the analogous final states in the large 
samples of vetoed J/ipK^ events. Events with misidentified muons escaping the D mass veto are accounted for using 
K^h ± fj, T events weighted by the per-particle probability, determined from uncorrelated PID control samples, for 
h to be misidentified. We also account for small contributions from misreconstructed signal events and charmonium 
events that escape the charmonium mass vetos. We test the fit using the large samples of vetoed J/ipK^*' and 
ijj(2S)K W events, and find good agreement with accepted values for branching fractions to the individual final states 
used here and the small charmonium isospin asymmetries 17| . 

We consider systematic uncertainties associated with reconstruction efficiencies; hadronic background parameteri- 
zation in di-muon final states; peaking background contributions obtained from simulated events; and possible isospin 
asymmetries in the background pdfs. We quantify efficiency-related systematics using the vetoed J/ipK^*' samples, 
including charged track, tt°, and reconstruction, PID, NN selection, and the AE and K* mass selections. Nearly 
all systematic effects largely cancel in the Aj ratio, and the final systematic uncertainties are small compared to the 
statistical ones. 

Table U shows the Af ( ' results. We find no significant isospin asymmetries in the high q 2 region, but we find 
evidence for large negative asymmetries in the low region. We calculate the statistical significance with which a null 
isospin asymmetry hypothesis is rejected using the change in log likelihood \/2A ln£ between the nominal fit to the 
data and a fit with Af ( ' = fixed. Figure [2] shows the likelihood curves obtained from the K£ + £~ and K*£ + £~ 
fits. The parabolic nature of the curves in the Af" > -1 region demonstrates the essentially Gaussian nature of 
our fit results in the physical region, and the right-side axis of Figure shows purely statistical significances based 
on Gaussian coverage. Incorporating the relatively small systematic uncertainties as a scaling factor on the change 
in log likelihood, the significance in the low q 2 region that Af is different from zero is 3.2cr for K£ + £~ and 2.7a for 
K*£ + £~ . We have verified these confidence intervals by performing fits to ensembles of simulated datasets generated 
with Af ' = fixed, and we find frequentist coverage consistent with the Aln£ calculations. The highly negative 
Af values for both K£ + £~ and K*£ + £~ at low q 2 suggest that this asymmetry may be insensitive to the hadronic 
final state, and so we sum the likelihood curves as shown in Figure M and obtain Af = -0.64±[J ^±0.03. Including 
systematics, this is a 3.9a difference from a null Af { ' hypothesis. 

For the angular analysis, because of the relatively small number of expected K*£ + £~ signal candidates in each 
q 2 region, a simultaneous fit over toes, cos 9k and cos^ is not possible, and an iterated fitting procedure is used. 
The first step is an toes fit, with the same underlying components as above, in each q 2 region for total signal 
(Ns) and combinatoric background (Nb) yields for the combination of all K*£ + £~ final states. The second fit is to 
co&Ok for events with toes > 5.27GeV/c 2 , where the only free parameter is Fl, and normalizations for signal and 
combinatorial background events are taken from the initial toes fit- The background normalization is obtained by 
integrating, for toes > 5.27GeV/c 2 , the ARGUS shape resulting from the toes fit- We model the cos 8k shape of the 
combinatorial background using e + e~ and events, as well as lepton-flavor violating e + /i _ and /i + e~ events, in 

the 5.20 < toes < 5.27GeV/c 2 sideband. Simulated events are used to account for the small remaining background 
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Table I: ' results. Errors are statistical and systematic, respectively. 
Mode low q 2 high q 2 



B^K^i+fi- -0.91±^ 2 ±0.18 0.39+°-^ ± 0.04 

•'-•< -i.4ilg:S = 

B -* Kl+i~ — 1.43±g;|| ± 0.05 0.28±°;^ ± 0.03 

B-*K*n+n- -0.26+°;34 ±0.05 -0.08t° 37 ■ 

B^K*e+e~ -0.66±%{% ± 0.02 0.32±°^| ± 0.03 

B->ir*l+^- -0.56l°;i5 ±0.03 nla+0 - 36 




contributions. The signal distribution given in Eq. 2 is folded with a model for the signal acceptance as a function 
of cos 9k obtained from simulated signal events. 

The final fit is to cos 6^, again for events with toes > 5.27GeV/c 2 , where the only free parameter is Afb- This 
fit requires the value of Fl as an input, which is taken from the second fit, as are the normalizations for signal and 
combinatorial background. The cos 9i shape of the combinatorial background is obtained using the same sideband 
samples as for the cos 9k fit. Correlated leptons coming from B -> £>M&/, D -> KW>Iv give rise to a peak in the 
combinatorial background at cos9i > 0.7 which varies as a function of toes, and we consider this variation in our 
study of systematic uncertainties. As for cos0£, the signal distribution given in Eq. 3 is folded with a model of signal 
acceptance as a function of cos 9i taken from simulated events. We again test our fits using the large sample of 
vetoed charmonium events, where Afb is expected to be zero, and branching fractions and the K* polarization are 
well-known 19|, |20(, and obtain results consistent with accepted values for each of the six individual K*£ + £~ final 



states, as well as their combination. We further test our methodology by performing the Toes and cos9i fits on a 
B + — ► K + £ + £~ sample, where Afb ~ is expected in most new physics models as well as the SM, and find values 
consistent with a null result in both low and high q 2 regions. 

We consider systematic effects from several sources. Uncertainties in yields due to variations in the ARGUS shape 
in the toes fits are propagated into both angular fits, and uncertainties on the fit Fl values are propagated into the 
Afb fits. Combinatorial background angular shapes are varied by dividing the sideband sample into two disjoint 
regions in Toes • We vary the signal model using simulated events generated with different form factors 12J, |21| and a 
wide range of values of CV, Cg and Cio- We perform fits to large ensembles of datasets for each generator variation to 
determine an average absolute fit bias. Finally, we constrain the modeling of mis-reconstructed signal events from the 
fits to the charmonium samples, where it is the largest systematic uncertainty. As with Af , the final systematic 
uncertainties are small compared to the statistical ones. 
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The results for Fl and Afb are shown in Fig. [TJ In the low q 2 region, where we expect Afb = —0.03 ± 0.01 [22J 



and F L = 0.63 ± 0.03 [10( from the SM, we measure A FB = 0.24±g;£§ ± 0.05 and F L = 0.35 ± 0.16 ± 0.04, where the 



first error is statistical and the second is systematic. In the high q region, the SM expectation is Afb ~ 0.26 and 

^;i±0.07andF i = 0.7ll°|°: 



F L ~ 0.40, and we measure Afb = 0.76^;^ ±0.07 and F L = 0.7lt°;22±0.04, with a signal yield of 36.6±9.6 events. 



Theoretical uncertainties on the expected SM Fl and Afb values in the high q 2 region are difficult to characterize, 
and the quoted values are obtained from our signal event generator 13, 21 1 . 

The magnitude of possible contributions from new physics to Cio can be constrained by a positive sign of Afb 
at high q 2 . By combining Afb with inclusive branching fraction results, an upper bound of |C U ^ P | <^ 7 can be 
obtained, improving on an upper bound derived solely from branching fraction results of |C^ P | <~ 10 j^ . The Afb 
results additionally exclude a wrong-sign C9C10 from purely right-handed weak currents at more than 3 standard 
deviations significance. The low q 2 AFB result suggests that a zero-crossing point, a distinctive SM feature, may 
not be present. Some angular asymmetries in the K* system include the longitudinal polarization as a component 
and are more sensitive to possible new physics contributions than measurements of Fl alone, as predicted values 
for these asymmetries in the SM and various new physics models can be calculated with relatively small theory 
uncertainties jSj ]. However, such asymmetries could not be considered here as they require analysis of a dataset 
substantially greater than currently available. Our results are consistent with measurements by Belle 25J, and replace 
the earlier BaBAR results in which only a lower limit on Afb was set in the low q 2 region 
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